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CALCULATION  OF  GROWTH  STRESS  IN  Si02  SCALES  FORMED  BY  OXIDATION  OF 
SiC  FIBERS 


R.  S.  Hay 

Air  Force  Research  Laboratory 

Materials  and  Manufacturing  Directorate,  WPAFB,  OH 
ABSTRACT 

A  numerical  method  to  calculate  growth  stress  in  Si02  scales  formed  during  SiC  fiber  oxidation 
is  described.  Calculations  were  done  for  SiC  fibers  between  700°  and  1300°C  using  previously 
measured  Deal-Grove  parameters  for  oxidation  kinetics  and  an  Eyring  viscoplastic  model  for  Si02 
scale  viscosity.  Initial  compressive  stresses  in  Si02  of  ~25  GPa  from  the  2.2x  oxidation  volume 
expansion  are  rapidly  relaxed  to  lower  levels  by  flow  of  silica  with  a  shear  stress-dependent  viscosity. 
At  700°  -  900°C,  axial  and  hoop  stress  at  the  GPa  level  persist.  Radial  expansion  of  the  outer  scale 
causes  hoop  stress  to  become  tensile;  axial  stress  becomes  tensile  by  the  Poisson  effect.  Tensile  hoop 
stresses  can  be  >2  GPa  for  thick  scales  fonned  at  <1000°C.  Effects  of  different  fiber  radii  on  growth 
stresses  are  examined.  Limitations  of  the  method  and  analytical  approximations  are  discussed. 

INTRODUCTION 

Constraint  of  the  2.2x  volume  expansion  during  oxidation  of  SiC  to  Si02  generates  very  large 
growth  stresses.  Microstructural  evidence  for  these  stresses  exists  for  crystalline  scales  on  SiC  fibers. 
High  dislocation  densities  in  crystalline  Si  02|  near  the  SiC-Si02  interface  suggest  high  shear  stresses 
exist  during  growth  of  new  crystalline  scale?  Axial  cracks  fonn  in  the  outer  scale  from  tensile  hoop 
growth  stress.1'3  Tensile  growth  stress  in  the  scale  may  decrease  fiber  strength  by  surface  nucleated 
fracture.4  Residual  stresses  affect  SiC  fiber  strength  and  are  important  for  SiC-SiC  ceramic  matrix 
composites  (CMCs)  for  high  temperature  structural  applications.5'6 

Growth  stress  during  silicon  oxidation,  which  has  a  volume  expansion  similar  to  that  for  SiC 
oxidation,  has  been  extensively  modeled.7"9  Most  recent  models  recognize  that  flow  at  high  stress  is 
non-Newtonian  (viscoelastic)  and  use  the  Eyring  model  for  shear-stress  dependent  viscosity.10'14 
Radial  compressive  and  tensile  hoop  growth  stresses  are  predicted  for  oxidation  of  silicon  fibers?2'15'16 
Axial  stresses  are  generally  ignored.10'1115'17  For  structural  fibers  this  is  an  important  omission, 
because  axial  stress  has  the  most  significant  effect  on  fiber  strength.  Applied  tensile  stress  increases 
oxidation  rates  of  silicon,18'20  and  recently  this  has  also  been  demonstrated  for  SiC  fibers.21 

A  method  to  calculate  the  radial,  axial,  and  hoop  growth  stress  components  anywhere  in  a  Si02 
scale  formed  by  oxidation  of  SiC  fibers  is  presented.  The  method  involves  discretization  of  the  scale 
into  layers  and  separate  calculation  of  stresses  in  each  layer.  Calculations  are  done  at  temperatures 
from  700°  to  1300°C  for  amorphous  scales.  Growth  stresses  for  SiC  fibers  with  different  radii  are 
examined.  Assumptions  and  limitations  of  the  method  are  discussed.  Complementary  results,  along 
with  a  more  thorough  description  of  the  background  information  and  method  are  given  elsewhere.22 

METHOD 

General 

A  schematic  of  the  volume  expansion  and  attendant  stresses  during  SiC  fiber  oxidation  with 
discretization  of  Si02  scale  into  annular  layers  is  in  figure  1.  There  are  two  sources  of  stress.  The  first 
is  the  initial  elastic  constraint  of  the  2.2x  volume  expansion  accompanying  oxidation  of  SiC  to  Si02. 
This  expansion  is  shown  in  stress-free  (dilational)  and  constrained  states  (Fig.  1).  The  second  is  the 
circumferential  expansion  of  old  scale  as  it  is  radially  displaced  outward  by  fonnation  of  new  scale 
(Fig.  1).  This  creates  tensile  hoop  stress  (ae)  in  the  outer  scale,4  23  and  tensile  axial  stress  (crz)  by  the 
Poisson  effect. 

The  large  volume  expansion  raises  concerns  for  growth  stress  modeling.  Elastic  constraint  of 
the  2.2x  volume  expansion  causes  ~25  GPa  compressive  stress  in  Si02  for  SiC  oxidation,  which  is 
much  larger  than  stresses  for  which  linear  elasticity  is  valid?  However,  high  shear  stresses  relax  very 
rapidly  to  values  appropriate  for  linear  elasticity  if  Si02  viscosity  is  shear-stress  dependent,  as  in  the 
Eyring  model,10'11’  '25  and  discretization  to  small  increments  can  confine  incremental  displacements 
and  stress  differences  to  values  appropriate  for  linear  theory.13’26'27 
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General  models  of  oxidation  growth  stress 
are  coupled  diffusion-reaction  and  fluid- 
mechanical  problems.26  The  oxygen  diffusion  rate 
drives  the  rate  that  stress  in  the  oxide  is  generated. 
Stress,  in  turn,  affects  the  oxygen  diffusion  rates. 
The  measured  oxidation  rates  for  SiC,  whether  they 
are  flat  plates,  fibers,  or  particles,  inherently 
include  the  effect  of  growth  stress;  stress-free 
oxidation  rates  cannot  be  measured.  The  effect  of 
stress  on  oxidation  rate  will  therefore  not  initially 
be  considered  when  modeling  growth  stress. 

The  following  assumptions  are  used  in 
modeling  the  growth  stress: 

1 .  Oxidation  volume  expansion  is 
dilational. 

2.  Stresses  resulting  from  constraint  of 
oxidation  expansion  are  relaxed  by  flow 
of  Si02  with  shear  stress-dependent 
viscosity. 

3.  Discretization  of  oxidation  to  small 
increments  allows  use  of  linear 
elasticity. 

4.  Growth  stress  effects  on  oxidation 
kinetics  are  not  considered. 

5.  Stress  relaxation  in  the  SiC  fiber  is 
negligible. 


Fig.  1.  Schematic  diagrams  of  oxidation  of  SiC 
to  SiOi-  The  original  SiC  surface  is  a  blue  line 
(c0,  bo),  midway  through  the  scale  thickness. 


Fiber  Oxidation  Kinetics 

Fiber  oxidation  kinetics  will  not  deviate  significantly  from  flat-plate  geometry  kinetics  until  the 
oxidation  product  for  12  pm  diameter  fibers  is  several  microns  thick.1,  5,28  SiC>2  thickness  (w)  (Fig.  1) 
therefore  obeys  Deal-Grove  kinetics  for  flat-plate  geometry:1,29 


4Bt 


A2 


1 


[1] 


The  parabolic  and  linear  rate  constants  are  B  and  B/A,  respectively.  B  and  A  obey  the  usual  Arrhenius 
relationships: 

-Qa 

A  =  A0e  RT  [2] 

-Qb 

B  =  B0e  RT  [3] 

where  T  is  absolute  temperature,  R  is  the  gas  constant,  QA  and  Qb  are  activation  energies  and  A0  and 
B0  are  pre-exponential  factors.  The  SiC  radius  (b)  after  oxidation  is  (Fig.  1): 


b  -  w2(Q2  —  Q)  +  b“5  —  Gw  [4] 

where  b0  is  the  original  fiber  radius  and  Q  is  the  ratio  of  SiC/SiCF  molar  volume  ratio.  The  outer  radius 
of  the  SiC>2  scale  (c)  is  (Fig.  1): 


c  =  b+w 


[5] 
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Elastic  Growth  Stress 

The  elastic  stresses  and  strains  are  determined  by  modification  of  a  method  used  for 
sequentially  deposited  coatings.  For  discretized  fiber  oxidation,  the  “deposition”  sequence  is 
reversed;  the  last  oxide  increment  fonns  at  the  SiC-SiO?  interface,  and  the  oldest  is  at  the  surface.  The 
system  is  divided  into  the  unoxidized  SiC  fiber,  the  SiC>2  added  from  time  t(i  -1)  to  t(i)  at  the  SiC-SiCF 
interface,  and  the  outer  SiCb  increment  added  fromct(0)  to  t(i-l)  (Fig.  1).  The  axial,  radial,  and  hoop 
stresses  in  the  SiC  fiber  at  time  t(i)  are  azSlC(i),  arSlCi[)  <  and  aeSlC(i)  respectively,  and  in  the  two  SiO 
increments  are  az  (i),  ar  (i)  ,  ere  (i)  and  az  (i-l),  ar  (i-l)  ,  ae  (i-l),  respectively.  The 
increment  in  SiCF  thickness  formed  from  t=i- 1  to  t=i  is  w(i)-w(i-l).  The  effect  of  the  i  Si02 
increment  on  growth  stresses  is  found  using  strain  compatibility  equations,  where  the  dilational  strain 


and 


[6] 


aziC(0  “ 


-Esicf 


Esicb2(i)  +  Esio2  ((b(i)  +  w(i)  -  w(i  -  l))2  -  b2(i)) 


TT 


01U9 

°z  (0 


CTsic(i)  =  apiC(i  —  i)  —  Pi, 

-SiC  _  _SiC/':  — 

CTe  (0  =  CTe  0  -  1)  -  Pi 
— EsioJ 


ti  ^ESicb(i)2  +  EsiC)2  ((b(i)  +  w(i))2  -  b(i)2)j 
Si02  s-\ 

(0  =  “Pi 


Si02,.,  Pj(b(0  +  w(i)) 

oe  20)  = 


Si02  /•  \ 

az  2(i  —  1)  = 


w(i)  —  w(i  —  1) 
f 


n  (b(i  —  1)  +  w(i  —  1)  —  w(i))2  —  b2(i  —  l)j 


af°2(i-l)  =  0 


Si02  /■  .  y* 

Gfi  (1  -  1)  = 


2pisb(i  -  l)2  -  Pi  ((b(i  -  1)  +  w(i  -  D)2  +  b(i  - 1)2) 
(b(i  —  1)  +  w(i  —  l))2  —  b(i  —  l)2 


+  CTziC(i  —  1)  [7] 

[8] 

[9] 

[10] 

[11] 

[12] 

[13] 

[14] 

[15] 


where  f  is  the  axial  force,  pjS  is  the  pressure  across  the  SiCF-SiC  interface,  pi  is  the  pressure  across  the 
interface  between  the  ith  and  the  i-l  h  Si02  increments,  ESic  and  ESi02  are  Young’s  modulus  of  the  SiC 
fiber  and  the  Si02  scale,  respectively,  and  vSiC  and  vSi02  are  Poison’s  ratio  for  the  SiC  fiber  and  the  Si02 
scale.  Stresses  in  older  increments  (j  =  i-2  to  j  =  0)  are  updated  with  the  stress  values  in  [13  -  15]: 

Relaxation  of  Elastic  Stress 

The  relaxation  of  the  elastic  stresses  for  all  increments  (j=l  to  i)  in  time  increment  At=t(i)-t(i-l) 
are  calculated  next.  The  Eyring  model  for  shear-stress  (x)  dependence  of  glass  viscosity  (r|)  is  used  for 
Si02:U"12’16’31"34 
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h  =  ho 


[16] 


TV, 


2kT 


Sinh  (T  Vc/ 2rt)  110  Sinh(T /tc) 


Vt 


Vc  is  the  activation  volume  for  plasticity  in  SiC^.1 1-12,34  k  is  Boltzmann’s  constant,  tc  is  the  critical 
shear  stress  above  which  plasticity  is  significant  (typically  — 100  MPa),  and  ri0  is  the  stress-free  SiCfi 
viscosity:35'36 

Pc  =  Co  ei^  [17] 


where  C„  and  Q  are  the  pre-exponential  factor  and  activation  energy  for  stress-free  viscosity,  respectively.35'36 
Shear  stress  relaxation  obeys  a  Maxwell  viscoelastic  model:1134,37 

^=_G7tl(T)'  X[t(i  -  1)]  =  T0  [18] 

where  G  is  the  SiC>2  shear  modulus,  and  the  initial  elastic  shear  stress  at  t(i  -1)  is  x0.  The  relaxation  of 
T0  to  a  new  value  (x)  in  time  increment  At  for  all  the  increments  (j=l  to  i)  is  determined  by  substitution 
of  [16-17]  in  [18]  and  solving  differential  equation  [18]  for  x: 


4kT 


T 


V 


-Coth" 


Gt 

e^o 


L\ 


Tanh 


VcTo(j) 

4kT 


=  2xcCoth  1 


Gt 

e^o 


Tanh 


x0(j) 


2t, 


[19] 


x0  is  determined  from  the  principal  stresses  for  all  the  increments  (j=l  to  i)  by  the  usual  method: 


38 


X0  (j)  =  [|  Sij  (j)Sij  (j)] 


/ 2 


[20] 


Sij(j)  are  the  components  of  the  deviatoric  stress  tensor  of  the  jth  annular  element  at  time  t(i): 


Sij  (j)  =  Oij  (j)  -  -Sjj  akk(j) 


[21] 


Relaxation  of  ae(j)and  cz( j)  is  proportional  to  x(j)/x0(j)  and  to  their  difference  with  ar(j),  which 
is  a  boundary  condition,  being  zero  at  the  SiC>2  surface  and  near  zero  elsewhere.  The  relaxed  values  of 
CTe(j)and  az(j)  (a0(j)'  and  <rz(j)')  are  determined  by  solution  of: 


ae(jy  0 

T(i) 

CTe(j)  -  ar(j)  0 

_i_ 

Gr(j)  0 

0  CTzGT. 

T0(j) 

0  nz(j)-Gr(j). 

i 

o  CTr(j). 

[22] 


Radial  Displacement  and  Hoop  Stress  Generation 

Relaxation  expands  the  SiC>2  scale  radially.  The  individual  radial  displacement  of  the  jth 
increment  (ur)  is: 


where 


ur(j)  =  - 

fisiC 
sii02(j)  = 


(i  +  Er02o)+4l02o)+4102a)4102o)) 

-t_ki02Q)  -  usl„  (trl10^)  +  °?02Q))1 

bSi02  L  v  /J 


[23] 

[24] 
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The  total  radial  displacement  of  the  jth  increment  is  the  sum  of  the  displacements  of  younger 
increments.  This  adds  hoop  strain  (eeSl0: *)  to  outer  layers  as  they  are  forced  to  a  larger  circumference 
(Fig.  1).  The  new  hoop  strain  in  each  increment  is: 


4i02ay  = 


i 

^  Ur(j) 


b(j  ~  1)  —  b(j) 
b(j) 


[26] 


Recalculation  of  Elastic  Stress  in  the  Scale  and  SiC  Fiber  after  Radial  Displacement 

The  stresses  in  each  S1O2  increment  are  recalculated  for  £eSl02(j)'  by  solving  the  three  strain 
compatibility  equations  for  the  three  principal  stresses.  Revised  axial  stress  czSlC  is  computed  from  the 
force  exerted  by  the  SiCL  scale,  which  is  the  sum  of  the  axial  stress  in  each  SiC>2  increment  x  area  of 
that  increment: 


1  x 


'  (j)  [2c(i)  (w(j)  -  w(j  -  1)  +  w(j)2  -  w(j  +  l)2] 


The  revised  radial  and  hoop  stress  in  the  SiC  fiber  are  calculated  by  determining  the  net  pressure  (pn) 
from  the  sum  of  the  pressures  in  each  annular  increment: 


1  x 

'(i)  =  a0lC(i)  =  -Pn  =  - 


These  revised  stresses  are  added  to  the  next  increment  i  as  the  program  loops  back  to  equations  [1  -28]. 
For  SiC  oxidation  when  w(i)  «  b(i),  the  stress  in  SiC  is  much  smaller  than  that  in  SiC>2. 


RESULTS  AND  DISCUSSION 
General 

Growth  stress  calculations  were  done  using  [1-28]  in  a  Mathematica™  program  using  scale 
discretization  to  500  layers  (i=500).  Calculations  were  done  for  oxidation  of  Hi-Nicalon™-S  SiC  fiber 
for  amorphous  scales  of  10,  100,  300,  1000,  and  3000  mn  thickness  (w)  at  700,  800,  900,  1000,  1100, 
1200,  and  1300°C.  The  Deal-Grove  oxidation  kinetics  for  this  fiber  have  been  reported  for  dry  air 
between  700  and  1300°C,  with  A0  =  6.5  x  10'4  m,  B0  =  1.2  x  10'8  nr/s,  QA  =111  kJ/mol,  and  QB  =  249 
kJ/mol.1,39  The  initial  fiber  radius  (bo)  is  6.1  pm.  The  molar  volumes  for  SiC  (Clsic)  and  amorphous 
Si02  (0Sio2)  are  27.34  cm3  and  12.46  cm3,  respectively.  The  Young’s  modulus  (E)  and  Poisson’s  ratio 
(v)  values  used  for  SiC  and  SiCL  were  Esic  =  400  GPa,  vSic  =  0.157,  ESio2  =  70  GPa,  and  vSi02  =  0.17. 
The  shear  modulus  for  silica  (G)  is  34  GPa,  and  has  only  weak  temperature  dependence.40  The  pre¬ 
exponential  factor  C0  and  activation  energy  Q  for  stress-free  SiC>2  viscosity  are  3.8  x  10’13  Pa-s  and 
712  kJ/mol,  respectively.35'36  The  activation  volume  for  plasticity  in  SiCL,  Vr,  decreases  with 
temperature,  and  has  been  inferred  to  have  values  ranging  from  1.2  x  1 0‘”  m  to  3  x  1 0"  nr.  '  ’  Vc 
corresponds  to  a  critical  shear  stress  xc  0 f  about  100  MPa,  which  is  roughly  consistent  with  experiment 
at  500  to  1400°C. 32,41  The  variation  in  growth  stress  with  change  in  fiber  radius  was  examined  by 
calculations  using  b0  =  2  pm  and  “flat-plate”  bo  A 00  (1  km).  Calculations  for  ae,  az,  a,,  and  x  at  T  of 
700°  -  1300°C  for  w  =  300  mn  and  b0  =  2  pm,  6  pm,  and  “flat -plates”  are  shown  ^in  figure  2. 
Calculations  for  other  scale  thicknesses  and  for  crystallized  scales  are  reported  elsewhere.22 

The  accuracy  of  growth  stress  calculations  is  limited  by  silica  viscosity  accuracy.  Silica  scales 
formed  during  SiC  oxidation  incorporate  carbon.42’45  Network  carbon  in  amorphous  SiCL  stiffens  the 
network  structure,  making  it  more  viscous  and  less  permeable  to  O2. 46-48  However,  SiCF  viscosity  may 
also  be  reduced  by  incorporation  of  impurities  in  Hi-Nicalon™-S  fiber.39 
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Flat  Plate 


r=  6  iim  r  =  2  \im 


Fig.  2.  Growth  stress  calculations  for  erg,  crz,  cr,-,  and  tat  T  of  700°,  800°,  900°,  1000°,  1100°, 

1200°,  and  1300°C  for  a  2  /um  radius  Jiber,  6  /um  radius  fiber,  and  a  flat-plate  for  a  300  nm 
amorphous  SiO 2  scale  on  Hi-NicaIon1M-S  SiC  fiber. 

General  features  for  hoop  stress  (ae)  and  axial  stress  (az)  in  amorphous  SiC>2  scale  are  evident 
in  figure  2  for  300  nm,  thick  scales  on  2  pm  and  6  pm  SiC  fibers  and  “flat-plate”  SiC.  The  continuous 
changes  of  ae  and  a,  with  change  in  bo  at  700°C  to  1300°C  throughout  the  scale  are  evident. 
Compressive  elastic  stress  of  — 25  GPa  for  ae  and  az  is  rapidly  relaxed  by  the  Eyring  model  [16]  for 
SiC>2  glass  viscosity  to  — 4  GPa  at  700°C  and  — 1  GPa  at  1200°C  just  1  nm  away  from  the  SiC-SiCF 
interface.  At  T  >1200°C  glass  viscosities  are  so  low  that  stress  levels  are  insignificant  more  than  20  nm 
from  the  SiC-SiC>2  interface.  For  flat  plates  ae  =  az,  and  stress  is  always  compressive.  az  is  driven 
towards  tensile  values  by  the  Poisson  effect  from  the  tensile  qe  caused  by  outward  radial  expansion. 
Note  that  a  300  nm  scale  takes  ~108  s  to  form  at  700°C  and  ~10  s  to  form  at  1000°C,  so  the  latter  is  of 
more  practical  interest.  Calculations  for  crystalline  SiC>2  scale  are  reported  elsewhere.  '  The  viscosities 
of  crystalline  scale  are  much  higher  than  those  for  amorphous  scale,  so  growth  stresses  in  crystalline 
scale  are  much  larger  than  amorphous  scale. 

For  a  2  pm  radius  fiber,  tensile  ae  develops  at  the  scale  surface  at  all  temperatures,  and  reaches 
3  GPa  at  700°C.  Tensile  az  does  not  develop  until  T  >  800°C  and  reaches  300  MPa  at  1000°C.  For  a  6 
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pin  radius  fiber,  tensile  Gq  does  not  develop  until  T  >  850°C  and  reaches  a  maximum  of  about  600 
MPa  at  1000°C;  tensile  az  does  not  develop  until  1100°C  and  reaches  100  MPa  at  that  temperature.  For 
thicker  scales,  larger  fractions  of  the  scale  are  under  tensile  stress.22 

For  a  2  pm  radius  fiber,  compressive  ar  up  to  -200  MPa  is  present  in  the  center  of  the  scale  at 
800°  -  900°C.  This  is  a  consequence  of  tensile  ae  at  the  scale  surface.  ar  decreases  towards  the  SiC- 
SiC>2  interface  from  the  high  compressive  ae  near  that  interface.  For  a  6  pm  radius  fiber  ar  is  tensile  at 
700°  -  900°C  because  ae  is  compressive  at  the  scale  surface  and  throughout  the  scale  thickness,  and 
only  reaches  small  compressive  values  at  T  >  900°C  when  significant  amounts  of  the  scale  are  in 
tensile. 

x  of  ~12  GPa  is  rapidly  relaxed  to  ~  2  GPa  at  700°C  and  500  MPa  at  1 100°C  just  10  nm  away 
from  the  SiC-SiCF  interface  (Fig.  2).  It  continues  to  decrease  as  compressive  az  and  ae  decreases  with 
distance  from  the  SiC-SiCF  interface  until  Gq  becomes  tensile,  at  which  point  x  increases  towards  the 
SiC>2  surface,  reaching  surface  values  >1.5  GPa  at  700°C  for  a  2  pm  radius  fiber  and  >  1  GPa  for  a  6 
pm  radius  fiber,  and  values  close  to  500  MPa  at  1000°C  for  both  fiber  radii.. 

Steady-State  Tensile  Stress 

For  the  300  nm  thick  scale  a  “steady-state”  develops  for  ae  tensile  stress  for  T  >  700°C  for  b0  = 
2  pm  and  T  >  1000°C  for  b0  =  6  pm.  For  thicker  scales  the  steady-state  region  develops  at  lower 
temperatures.22  For  x  <  xc  of  100  MPa  [16]  where  a  stress-free  viscosity  |r|0)  [17]  is  applicable,  an 
analytical  expression  for  the  steady-state  hoop  stress  [ae(ss)]  can  be  derived:" 


Ge  (ss)  = 


[29] 


For  x  >  100  MPa  the  stress  dependence  of  viscosity  is  significant,  and  r|  [16]  must  be  substituted  for 
x\o,  giving: 


ct@  (ss)  =  — 2xcCsch 


bwxcn 
■Bti0  (ft  -  1). 


[30] 


Comparison  of  predictions  for  steady-state  ae(ss)  tensile  stress  from  [30]  for  a  6  pm  radius  (b0)  fiber 
with  numerical  calculations  are  shown  to  be  very  close  in  another  publication.22 

By  insetting  the  Arrhenius  expressions  for  B  and  r|0  into  [29],  the  decrease  in  steady-state 
tensile  stress  with  increasing  temperature  is  a  consequence  of  Q  >  Qt,: 


Ge  (ss)  = 


Q-Qb  / 1 
2B0C0e  RT 

wb 


0 


[31] 


If  materials  exist  for  which  Q  <  Qb,  an  increase  in  tensile  growth  stress  with  increase  in  oxidation 
temperature  is  expected. 

Summary  and  Conclusions 

A  method  to  calculate  the  axial,  hoop,  and  radial  growth  stresses  in  SiCF  scales  generated  by 
the  2.2x  volume  expansion  during  SiC  fiber  oxidation  was  developed.  The  method  assumes  the  initial 
oxidation  volume  expansion  is  equal  in  all  directions  (dilatational)  and  that  the  stresses  resulting  from 
constraint  of  that  expansion  are  relaxed  radially  with  an  Eyring  stress-dependent  SiC>2  viscosity, 
although  other  appropriate  viscosity  models  can  be  substituted.  The  method  can  be  equally  well 
applied  to  fibers  of  silicon  or  other  materials.  High  compressive  hoop  and  axial  stresses  of  ~25  GPa  are 
very  quickly  relaxed  to  much  lower  values  at  all  temperatures.  Radial  expansion  creates  tensile  hoop 
stress  in  the  outer  scale.  Tensile  hoop  stress  eventually  drives  axial  stress  to  a  tensile  state  by  the 
Poisson  effect.  Tensile  hoop  and  axial  stress  can  reach  values  >  2  and  0.5  GPa,  respectively  for 
oxidation  for  long  times  at  700°  -  900°C  on  Hi-Nicalon™-S  fibers.  At  temperatures  greater  than 
1200°C  growth  stresses  are  quickly  relaxed  to  negligible  levels  by  viscous  flow  of  Si02.  The  accuracy 
of  the  growth  stress  calculation  method  is  likely  to  be  limited  by  knowledge  of  accurate  values  for 
amorphous  silica  viscosity.  Tensile  hoop  stresses  reach  steady-state  values  that  can  be  described  by 
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analytical  expressions.  The  decrease  in  tensile  hoop  growth  stress  with  increase  in  oxidation 

temperature  is  a  consequence  of  activation  energy  for  viscous  flow  >  activation  energy  for  oxidation. 
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